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ABSTRACT 



o 

O Aims. We identify a prominent absorption feature at 1 1 15 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret 
its astrophysical origin. 

C/j Methods. The characteristic hyperfine pattern of the H20^ ground-state rotational transition, and the lack of other known low-energy transitions in 
, ^ , this frequency range, identifies the feature as H2O* absorption against the dust continuum background and allows us to derive the velocity profile 
of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the 
I frequency of the transition and the conditions for its formation. 
^ Results. In DR21, the velocity distribution of H20^ matches that of the [Cii] line at 158 /jm and of OH cm-wave absorption, both stemming 
I from the hot and dense clump surfaces facing the H ii-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates 
the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H20^ column density of 
7.2 X 10'- cm-2 in NGC 6334, 2.3 x lO'"* cm"- in DR21, and 1.1 x 10'^^ cm^^ in Sgr B2. 
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1. Introduction detected. Section 3 summarises the spectroscopic data of the 

I ' J 1 1 . molecule. The observations and the line identification are de- 

^ ; Oxidaniumyl or oxoniumyl d Connelly et al.i |2005|), the reactive g^^j^ed in Sect. 4 and in Sect. 5 we discuss the physical proper- 

• ^ water cation, H20^, plays a crucial role in the chemical net- ties of the HiO^ absorption layer 
work describing the formation of oxyg en-bearing molecules iii 
H I UV irradiated parts of m olecular clouds jvan Dishoeck & BlacE 

! [T986HGerinetalil20Toh . It was identified at optical wavelengths 2 The SOUrceS 
in the tails of comets in the 1970's (Fehrenbach & Arpignvl 

ll973':'Herzb erg & Lewl ll974-, Wehinger et al., 1974), but its de- We observed three massive Galactic star-forming/H 11 regions 
tection in the general interstellar medium has proven elusive. with very different properties. The DR21 star-forming region is 
We report a detection of the ground-state rotational transi- embedded in a ridge of dense molecular material that obscures 
tion of HaO^ in some of the first spectra taken with the HIFI in- it at optical wavelengths. The embedded cluster drives a violent 
strument (de Graauw et al 2010) on board the Herschel Space bipolar outflow and creates bright photon-dominated (or photo- 
Observatory (Pilb ratt et all I2OIO ) during the performance ver- dissociation) regions (PDRs), visible as clunips of 8 ;um PAH 
ification campaign and early science observations. Section 2 emission in Spitzer IRAC maps (iMarston et al.Ll2004l) and show- 
briefly introduces the properties of the sources where HaO^ was ing up in emission lines from tracers of irradiated hot gas, such 

as HCO "'', high-7 CO, atomic and ionised carbon, and atomic 

* Herschel is an ESA space observatory with science instruments oxygen (|Lane et al., 1990; Jakob et al., 2007). The eastern, blue- 
provided by European-led Principal Investigator consortia and with im- shifted outflow expands in a blister-like fountain, while the west- 
portant participation from NASA. ern, red-shifted outflow is confined to a small cone. 
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Fig. 1. Energy level diagram of the lowest rotational levels of 
ortho-H20^ and its radiative transitions. The fine structure tran- 
sition frequencies are given in GHz. 



Table 1. Parameters of the hyperfine lines F' - F" in the ob- 
served 1 11 - Ooo, J - 3/2 - 1/2 ortho H20^ transition, including 
predicted frequencies, Einstein-A and optical depth at low tem- 
peratures 
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" Predictions based on IStrahan et all ( Il986l) and iMiirtz et alj ( Il998h . 
Nominal uncertainties are « 2 MHz but this is inconsistent with the 
discrepancy between the two predictions so that the actual uncertainty 
is unknown. 

from the matching DR21 OH pattern bv lGuilloteau et al.l j 19841) 

J'Tiow-r/A'HiO+'iv in lO"'* km s"' cm^ 

Table 2. Summary of the observational parameters 



The Sgr B2(M) and (N) cores are the most massive star- 
formation sites in our Galaxy. The line of sight, located in 
the plane of the Galaxy, passes through many spiral arm 
clouds and the extended envelope of Sgr B2 itself. The fore- 
groun d clouds display a very r ich molecular and atomic spec- 
trum (iPolehampton et aU |2007|) . although they often have very 
low densities and column densities, characteristic of diffuse 
or translucent clouds. The envelope of Sgr B2 itself includes 
hot, low density layers at both the ambient cloud vel ocity of 
64 km s and at km s ' (iCeccarelU et al.L l2002h . Many 
species detected along this line of sight have not been found 
elsewhere and the exact origin of the molecular features is of- 
ten a mbiguous because of the overlapping radial velocities (e.g., 
ICoinitoet^ll 120031) . 

NGC6334 is a nearby molecular cloud complex containing 
several concentrations of massive stars at various stages of evo- 
lution. The far-infrared source "f contains an embedded clus- 
ter of NIR sources (Tapia et al., 1996). Four compact mm con- 
tinuum sources are locate d near the geometric centre of the 
cluster (iHunteretal.Ll2007h . Although NGC6334I is not known 
to exhibit strong absorptio n lines, its OH absorption profiles 
(iBrooks & WhiteoakL I200TI) reveal two molecular clouds along 
this line of sight, one with velocities between - 1 5 and 2 km s ' , 
and the other near 6 km s ' . 



3. The H2O+ spectroscopy 

H20^ is a radical with a electronic ground state and bond 
lengths and angle slight l y larg er than H2O. Quantum-chemical 
calculations ( Weis et all 1198 91 yield a ground-state dipole mo- 
ment of 2.4 D. The B\ symmetry of the ground electronic state 
leads to a reversal of the ortho and para levels relative to water. 

The rotational spectrum was measured by laser magnetic res- 
onance (Strahan et al., 1986; Miirtz et al., 1998). Predictions of 
the Nk„k, = 111 - Ooo, -/ = 3/2 - 1/2 fine s tructure cornponen t 
near 1115 GHz using the new parameters bv lMiirtz et alJ 119981) 
are b etween 27.3 and 28 .5 MHz higher than those calculated 
from IStrahan et all (|I986'), even though both articles claim to 
have reproduced the experimental data to ~2 MHz. The reanal- 
ysisof equivalent measurements of SH^, by iBrown & Miilleii 
(i2009h . shows that this accuracy is in principle achievable. 
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However, the large centrifugal distortion in H20^ requires a 
large set of spectroscopic parameters to reproduce a compara- 
tively small set of data; this may cause problems in the zero-field 
extrapolation. Moreover, the frequencies of the tw o fine struc- 
ture levels of the In rotational state in Table V of iMiirtz et al.l 
(1998) agree precisely with those of the F' - J' , F" - J" hyper- 
fine transitions. This can only be achieved when the calculated 
frequencies are lower by 51.56 and 88.05 MHz, respectively, 
since the respective hyperfine component is the lowest in each 
case. Correcting the published frequencies of the J - 3/2 - 1/2 
fine structure component by 51.56 MHz improves the agree- 
ment with^trahan et al. (1986). The results are summarized in 
Table [H Alt ernatively, we could use the corrected frequencies of 
iMiirtzet alj ([r99&) and arrive at values that are lower by about 
23 MHz. This provides a rough estimate of the uncertainty in the 
predictions. An H20 ^ catalogue entry will be prepared for the 
CDMS (Mull er et all 12005) bv carefully scru tinizing the avail- 
able IR data summarised in lzheng et alJ(l2008h with x 150 MHz 
uncertainties. 



4. Observations of the 1115 GHz ground-state 
transition 

The H20^ line was detected in DR21 during performance ver- 
ification observations of the HIFI instrument, testing spectral 
scans in the HIFI band 4b. Later science observations of Sgr B2 
and NGC 6334 also confirmed the detection in these sources us- 
ing the identification and frequency assignment from DR21. The 
main parameters of the observations are summarised in Table |2] 
At 1 1 15 GHz, the Herschel beam has 21" HPBW. 
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Fig. 2. Fit of the hyperfine multiplet of the H2O+ 1115 GHz line 
in DR21. The bottom panel shows the 0.5 K absorption line 
superimposed on two different fit profiles, one based on a 3- 
component Gaussian (see text) and the other one usin g the OH 
6cm absorption spectrum from Guill oteau et alj (|1984|) . The top 
panel shows a breakdown of the fitted profile into its hyperfine 
constituents in the case of the 3-Gaussian profile. 



The identification with H20^ was straightforward in DR21 
because of the simple source velocity structure that cannot be 
confused with the well resolved, characteristic hyperfine struc- 
ture of the line. When fitting the line, one has to take into account 
that the line extinction begins to saturate, with a maximum opti- 
cal depth of 0.59 for DR21 and 1 .55 for Sgr B2 (see below). For 
DR21, we fitted the observed profile using an adjusted velocity 
profile with asymmetric wings. Because of the limited signal- 
to-noise ratio, the fit was performed manually by adding three 
Gaussian components of increasing width (see Fig.|2l). 

The resulting velocity distribution allows us to interpret the 
origin of the absorbing material by comparing with the veloc- 
ity distribution of other species observe d towards the same po- 
sition with comparable beam size (see Ossenkopf et al.L l2010t 



sition with com parable beam size (see UssenKopi et al.L UUlit 
iFalgarone et"!!! boiOtlvan der Tak et al.L I2OIOI) . Figure [3] shows 



that the peak HaO^ velocity of -1.7 km s~ is not seen in any 
other tracer. The intrinsic velocity of the DR21 molecular ridge 
is -3.0 km s which is matched by the line centres of the 
H'^CO^ 1-0, the CO 6-5, and the '^^CO 6-5 transitions. The 
higher excitation lines of '^CO, C"^0, H2O, and the [Cii] Hne 
exhibit a slightly blue-shifted peak velocity of about -5.0 km 
s The HaO^ profile exhibits a prominent, very broad blue 
wing. This is not present in any of the molecular emission lines. 
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Fig. 3. Comparison of the fitted H20^ velocity profile to other 
tracers observed in DR21 with similar beam size. The pro- 
files are normalised to a peak of unity and separated by multi- 
ples of 0.5 fr om bottom to top. The fit (bottom line) used the 
IStrahan et al.l ( 1986) based line frequency prediction, the profile 
at the top is shifted by -4.0 km s corresponding to a rest fre- 
quency lower by 15 MHz. 
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Fig. 4. Fit of the observed H20^1ine in Sgr B2. The dashed line 
visualises the velocity structure of the absorbers by plotting the 
strongest hyperfine component on a linear column density scale, 
i.e., without optical depth correction. 



but is found in the [C iil profile and t he OH absorption spectrum 
measured bv' Guilloteau et all (|1984|) towards the same position. 

To underline this good match, we have superimposed in 
Fig. [2] the absorption profile that would be obtained by simply 
performing the hyperfine superposition of the 6.030 GHz OH 
absorption profile. The match is as good as that achieved with 
the analytic profile and even reproduces the small excursions at 
1115.22 and 1115.27 GHz. This indicates that OH and H2O+ 
occur in the same region and under the same physical condi- 
tions. The displacement of the fitted profile relative to the [Cii] 
and OH profiles of about 4.0 km s ' is within the discrepancies 
between the different predictions of the line frequency. The as- 
tronomically determined line rest frequencies from comparison 
with the OH line fall 15 MHz below the predicted frequencies. 
As the line peak is very sharp, the accuracy of the frequency is 
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Fig. 5. Same as Fig.g] but for NGC6334. 



probably better than 2 MHz. Assuming a match with the [C ii] 
line instead, would provide a larger uncertainty of x6 MHz. 

The identification and the corrected frequencies are then 
used to analyse the line structures in Sgr B2 and NGC 6334 
(Figs.|4]and|5]l. In Sgr B2, we see absorption at both the velocity 
of its envelope and the velocities of many foreground clouds, al- 
most saturating the line. NGC6334 exhibits weak H20^ absorp- 
tion at -13 km s This deviates f rom the OH absorption pro - 
file towards the source measured bv lBrooks & Whiteoakl(l200lh . 
At velocities below -10 km s only some OH maser emission 
was found. This might indicate that the observed H20^ is not 
related to the foreground material, but to hot gas in the direct 
vicinity of the continuum s ources. Alternativel y, if we use the 
predicted frequencies from I Strahan et al.l (Il986l) in Table \T\ the 
H20^ absorption in NGC6334 is centred on -9 km s \ in rea- 
sonable agreement with the OH absorption at - 8.2 km/s mea- 
sured toward component fQ- At about -9 km s~' . iBeuther et"aLl 
(l2005h also observed CH3OH and NH3 absorption towards the 
H n region. 

5. Discussion and outlook 

That HiO^ shows up in absorption against the dust continuum 
implies that the excitation of the molecule must be cold er than 
the dust. As a reactive io n (see the discussion by Blaclj |2007t 
IStauber & Br uderer 2009' for C0+), H2O+ is not expected to be 
in thermal equilibrium at the kinetic temperature of the gas. Its 
excitation reflects either the chemical formation process or the 
radiative coupling with the environment. From a single absorp- 
tion line, one can only provide a lower limit to the H20^ column 
density, assuming a low excitation temperature where basically 
all H20^ resides in the ground state, which is applicable to tem- 
peratures well below the upper level energy of 53 K. 

Table[T]provides the integral over the optical depth of the hy- 
perfine components in the low temperature limit. For the overall 
J - 3/2-1/2 fine structure transition, we obtain a line integrated 
optical depth of / TdvlN^^o* = 4.70 10"'^ km s"' cm^ per 

A similar case is reported by ' Gerin et all ( |201(I ) for W31C. The 
source shows a complicated spectrum with multiple absorption compo- 
nents, but a closer co rrelation with other tracers is found when using the 
IStrahan et al.l l ll986h based frequencies. A recent detection of H20^in 
W3 IRS5 and AFGL259 1 by Benz et al. (i n prep.) seems to favour the 
frequency predictions bv lMiirtz et al.l ( 119981) . 



molecule, resulting in a lower limit to the H20^ ground-state col- 
umn densities of 7.2 x lO'- cm"^ for NGC 6334, 2.3 x lO'^ cm"^ 
forDR21,and 1.1 x 10'^ cm ^ for SgrB2. 

These values are lower limits not only because of to the low- 
temperature approximation, but also because they assume that 
the absorption occurs in front of the continuum source and not 
within the dusty cloud, where the line absorption is partially 
compensated by dust emission. There may also be additional 
amounts of H20^ in the para species that would not contribute 
to the 1115 GHz line. Altogether, the total H20^ column density 
could be much higher than the lower limits given here. 

The excellent correlation between the H20^ profile and the 
OH absorption profile in DR21 i ndicates that both s pecies occur 
in the same thin layer of hot gas ^ Jones et al.Lll994l) that directly 
faces the H 11 region at the blue-shifted blister outflow. There is 
no obvious coiTelation with the distributions of CO, H2O, or 
HCO^. For Sgr B2, we can clearly identify absorption in mul- 
tiple translucent foreground clouds. Their densities must be high 
enough to produce some molecular hydrogen, but low enough 
not to quickly destroy the H20^. For NGC6334, the gas compo- 
nent producing the H20^ absorption remains unidentified. 

With the identification of H20^ in the interstellar medium, 
we provide a first step to quantifying an important intermediate 
node in the oxygen chemical network, connecting OH^ in diff'use 
clouds and at cloud boundaries, through HjO^, with water in 
denser and cooler cloud parts. To obtain an estimate for the total 
H20^ abundance, we need to measure the excitation temperature 
of H20^. Observations of additional transitions of H20^, such as 
those at 742 GHz, are therefore essential. 
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